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quantum Monte Carlo simulation focused on the critical behevior near the quantum critical 
point. 
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1. Introduction 

In this article we discuss phenomena observed in sys- 
tems with a singlet ground state with triplet excitations 
separated by a finite spin gap. The focus of researchers' 
attention has been on the magnetic phase transition in- 
duced by a magnetic field and/or a pressure. While the 
field- (or pressure-) induced magnetic ordering were ob- 
served in various materials, the character of the phase 
transition, in particular the one at the quantum critical 
point (QCP) seems to be insensitive to the the nature of 
the ground state and to the origin of the gap. In one class 
of materials, localized magnetic spins are bound strongly 
in pairs to form singlet dimers, and the inter-dimer in- 
teraction is not strong enough to induce any long-range 
magnetic ordering. In these materials, the excited states 
consist of excitations of quasiparticles, recently often re- 
ferred to as tripletons or triplons. They originate from the 
triplet states at each site. The magnetic ordering in the 
dimer systems with a sufficiently strong magnetic field 
was suggested based on a mean-field theory. 1 ' 2 Later, two 
possibilities were suggested 3 depending on the relative 
strength of the repulsive interaction among the quasipar- 
ticles and their kinetic energy. If the kinetic term dom- 
inates, the system may undergo a condensation transi- 
tion as suggested by the mean-field theory 1,2 and later 
also by the Hartree-Fock (HF) theory with the bosonic 
representation. 4 On the other hand, when the repulsive 
force dominates, the ordered state is the super lattice 
where the quasiparticles are located periodically to form 
a lattice with a larger unit cell than that of the origi- 
nal lattice, which gives rise to magnetization plateaus in 
the magnetization-field curve. This scenario is realized in 
materials for which the effective dimer-dimer couplings 
vanish due to the cancellation among the exchange cou- 
plings and higher order interactions play essential roles. 5 
A well-known example is SrCu2(B03)2- 6 ~ 8 In the rest 
of the present article, however, we only discuss the first 
case, i.e., the case where the kinetic-term is dominant. 

2. Early Experiments and Mean-Field Theory 

A long lasting flow of publications on the field-induced 
magnetic ordering was initiated by a theoretical work, 1,2 
which was inspired by experiments 9, 10 on a quasi one- 
dimensional magnetic compound Cu(N03)2 • 2.5 H 2 0. 



The zero-field susceptibility measurement 11 had been 
suggesting that the system's main magnetic structure 
consisting of isolated dimers. Heat capacity measure- 
ments 12 and magnetization measurements 13 in the pres- 
ence of a magnetic field were also carried out. How- 
ever, no clear indication of a phase transition was ob- 
served. In retrospect, this was simply because the tem- 
perature was not low enough or the magnetic field was 
not strong enough to bring the system into an ordered 
state. Later, an unusual magnetocaloric behavior was ob- 
served, 9 which seemed to indicate some sort of phase 
transition. Subsequently, the existence of the transition 
was confirmed by an anomaly in the adiabatic magneti- 
zation curve 10 at 3-4 T below 1 K. 

In order to understand this phenomena, a model sys- 
tem that consists of dimers of spins was considered. 1,2 
The building block of the model is a pair of S = 1/2 
spins coupled to each other antiferromagnetically. When 
a magnetic field is applied to the system, the degeneracy 
among the triplet states is lifted and the lowest branch 
comes down to meet the singlet state. To be concrete, 
the spin Hamiltonian considered was 



Tt = 22 Jsg-Si2- 



(ij) a.,0 ia 



(1) 



Here J > is the antiferromagnetic coupling constant 
that binds a pair of spins, i and j the indices of dimers, 
and a and (3 the specifiers of individual spins in a dimer. 
The symbol Si a stands for an S = 1/2 spin operator. The 
dominant couplings are the intra-dimer couplings repre- 
sented by J. While the compound Cu(N0 3 ) 2 -2.5 H 2 
seemed to have a quasi-one-dimensional lattice struc- 
ture, 14 any specific structure was assumed in the mean- 
field theory because it is not relevant in the approxima- 
tion. In fact, as long as the lattice is three dimensional, 
its structure does not affect asymptotic behaviors near 
the criticality even in more elaborate treatments. There- 
fore, in most of the following discussion we only assume 
that the lattice is three dimensional. 

The above Hamiltonian can be further simplified, tak- 
ing into account the fact that we can neglect local high 
energy levels, when the inter-dimer couplings are suffi- 
ciently small and can be treated as a perturbation. Only 
two lower levels, the singlet state and one of the triplet 
states are retained. These two levels can be regarded as 
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the up and down states of an effective 5=1/2 spin. 
Then we have the following effective Hamiltonian: 1 

H = 4 (s?Sf + SfS] + l -S*S*\ Y, HiSf (2) 

(ij) V J i 

where is an operator representing the effective spin. 
The coupling and Hi can be expressed as simple linear 
combinations of J, H and J' ia jf 3 - 

We here consider the case where the effective Hamil- 
tonian is translation invariant and non-frustrated. The 
former is the case with most of materials unless impuri- 
ties affect the observations. The latter assumption does 
not strictly hold. However, in many important cases, the 
frustration is weak, at least not strong enough to alter the 
phase diagram structure or the critical behavior qualita- 
tively. With these assumptions, we obtain an 5 = 1/2 
XXZ model with easy plane anisotropy. The mean- 
field approximation was applied to this model, yielding 1 
Wmf = £(«) Jij (Sfm x + \Sfm z ) - H t m z . By solv- 
ing the self-consistent equation in terms of m x and m z , 
a symmetry-breaking solution was obtained, from which 
the presence of a field-induced magnetic ordering in real 
materials was suggested. 

However, in subsequent works it was found that the 
anomalous behavior observed in the adiabatic magne- 
tization curve 10 of Cu(N0 3 ) 2 ■ 2.5 H 2 should be at- 
tributed to the development of a short-range order along 
the chain, not a true long-range order. By using the exact 
solution 15, 16 to the 5=1/2 XY model in one dimension, 
it was shown 17 ' 18 that a purely one-dimensional model 
can produce the adiabatic magnetization curve similar to 
the one observed in the experiments. Later, this result 
was confirmed by a more extensive investigation. 19 An 
observation of the genuine long-range magnetic ordering 
had to wait till the NMR experiment, 20 which revealed 
that the emergence of the long range order takes place in 
two stages; first the development of a short range order 
along the chain at an intermediate temperature and then 
the three-dimensional long range order at a lower tem- 
perature. The heat capacity measurement 21 showed these 
characteristics more clearly in the form of two peaks; a 
broad heat capacity peak at a high temperature and a 
sharper one at a low temperature. The low-temperature 
one is identified as the transition discussed in the pre- 
ceding theoretical works 1,2 based on the mean-field the- 
ory. Subsequently, the magnetic structure in the ordered 
phase was determined through a neutron diffraction ex- 
periment 22 and an NMR experiment; 23 weakly coupled 
magnetic chains with each chain having alternating cou- 
pling constants. 

It should be pointed out here that the field-induced 
magnetic ordering itself is a rather generic phenomenon 
and not restricted to the dimer systems. For example, 
a theoretical result similar to the one mentioned above 
was obtained for an 5 = 1 (not necessarily dimer) spin 
system with an uniaxial spin anisotropy. 24 Due to the 
anisotropy, the single-site ground state is the S z = 
state whereas the S z = ±1 states are degenerate ex- 
cited states. The mechanism of the magnetic ordering 
induced by a magnetic field is quite analogous to that 



in the dimer system. A heat capacity measurement was 
carried out 25 on this type of material. The compound 
was Ni(C 5 H 5 NO) 6 (C104)2, in which each Ni 2+ ion car- 
ries an 5 = 1 spin. The spins are influenced by the uni- 
axial spin anisotropy of which the magnitude is greater 
than the spin-spin exchange couplings. A cusp-like peak 
was observed in the temperature dependence of the heat 
capacity, and the parabolic phase diagram in the H-T 
plane was drawn. There are several other reports on 
Ni based compounds. A magnetization measurement 26 
was performed for NiCi2-SC(NH2)2- Later, heat ca- 
pacity measurements on Ni(C5Hi4N 2 )2N 3 (C104) 27 and 
Ni(C 5 Hi4N 2 )2N 3 (PF 6 ) 28 clearly showed a lambda peak 
in the heat capacity. Compounds with more complicated 
spin anisotropics were discussed within the mean-field 
theory 29 and the spin-wave theory. 30 Recently, the rare- 
earth compound PrOs 4 Sbi 2 was studied through the 
heat capacity measurement, 31 and the neutron diffrac- 
tion. 32 The level crossing scheme of this material was 
theoretically clarified. 33 

The Cr-based compound Cs3Cr 2 Brg may belong to 
another category of materials that show a field-induced 
magnetic ordering. This is a dimer system in which spins 
are bound pairwise antiferromagnetically to form dimers. 
The Cr 3+ ions carry 5 = 3/2 spins. What makes this 
compound unique is the fact that the ions form frustrated 
triangular lattices stacked on top of each other. Due to 
the frustration, the magnetic ordering in this compound 
seems to be more complicated than a mere Neel state. A 
cusp in the temperature dependence of the susceptibility 
was observed 34, 35 in this compound although the nature 
of the ordered phase was not clarified. Corresponding 
to the four possible values of dimer angular momentum, 
5 = 0, 1, 2, 3, four plateaus were observed in the field de- 
pendence of magnetization at low temperatures. 36 Quite 
recently elastic and inelastic neutron scattering experi- 
ments 37 suggested that the gap does not close between 
two subsequent plateaus, in contrast to the other materi- 
als mentioned above. Theoretical investigations have not 
been done yet on this material. 

3. Recent Experiments and Bose Gas Represen- 
tation 

After a relatively dormant period of mid-90's, active 
researches of the field-induced ordering were restarted 
on copper based dimer compounds. The compound 
Cu2(C5Hi2N 2 )2CLi, in which Cu 2+ ions form weakly- 
coupled ladders, was investigated by various measure- 
ments; the magnetization, 38 the ac susceptibility, 39 the 
heat capacity 39 and the inelastic neutron scattering. 39 
The heat capacity measurement at higher magnetic field 
was performed later. 40 The sharp cusp was observed at 
the transition temperature and a parabolic curve of the 
critical temperature as a function of the magnetic field 
was drawn. 

Several other materials were also investigated. A 
heat capacity measurement showed a clear anomaly in 
(CH3) 2 CHNH 3 CuCl3. 41 An 5 = 1 dimer compound 
Ba 3 Mn 2 0s, in which Mn 5+ ions carry 5 = 1 spins and 
form singlet dimers, was studied 42 to reveal two step 
magnetization curve, corresponding to the three possi- 
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ble values of the dimer magnetization. 

Even more thoroughly investigated was the family of 
compounds XCUCI3 where X = K, Tl, NH4, etc. From 
susceptibility measurments 43 the broad peak tempera- 
ture T p was obtained for KCUCI3 and TICUCI3 at zero 
field. For these two compounds, the magnetization pro- 
cess was investigated at finite external field, 44 which 
yielded estimates of the critical magnetic field and the 
gap. An inelastic neutron scattering was carried out for 
KCuCl 3 45 and later for TlCuCl 3 , 46 which showed a three- 
line structure in the magnon excitation spectrum corre- 
sponding to the three triplet excitations. 45 By comparing 
the experimentally determined dispersion with the RPA 
approximation results, the exchange couplings were also 
estimated. 45 

The experimental evidences were accumulated and it 
was realized that, apart from the correct critical expo- 
nents for the phase transition, there are several addi- 
tional qualitative features that could not be accounted 
for by the mean field theory; 1 ' 2 (i) the upward convex 
curvature in the adiabatic magnetization curve, (ii) the 
upward convex curvature in the magnetization curve as 
a function of the temperature with fixed magnetic field, 
and (iii) the asymptotic behavior of the critical magnetic 
field as a function of the temperature that seems more 
like an algebraic dependence than an exponential depen- 
dence. 

The upward convex curvature in the adiabatic mag- 
netization had been known from the very early period 
of the study. 10 The temperature dependence of the mag- 
netization was studied for TICUCI3 47 and, at sufficiently 
strong fixed magnetic fields, revealed an upward convex 
curvature and a dip at a temperature slightly higher than 
the transition point. The critical magnetic field H C \(T) 
was estimated for various temperature and the fitting to 

H c i (T) — H c i (0) oc (3) 

yielded <j> — 1.7(1). 47 For the estimate of <fi, there are 
variations: = 2.2(1), 4 = 2.0(1), 48 = 2.3(1) 49 
and <\> = 1.67(7) 50 for TlCuCl 3 ; and = 2.6(2) 51 ' 52 for 
NiCi2-4SC(NH 2 )2- It is notable that all the estimates so 
far are higher than the mean-field prediction <fi = 3/2 as 
discussed below. 

The algebraic dependence (3) is in contrast to the 
mean-field theory which indicate an essential singular 
behavior near zero temperature: H c i(T) — H c i{0) oc 
exp(— Tq/T) with T being a constant. At least this in- 
consistency can be removed by introducing a mean-field 
theory (the HF approximation) based on the bosonic rep- 
resentation 4 of the Hamiltonian. More importantly per- 
haps, it provides a simple picture for the mechanism of 
the field-induced magnetic ordering; the BEC of triplet 
quasiparticles. This theory is equivalent to the spin wave 
theory with the M z — state taken as the vacuum 
state. A similar spin wave theory was considered in an 
early study. 17 It was found that the spin wave theory 
yields properties closer to the experimental observations 
than the mean-field theory in the spin representation. 
For example, the spin wave theory correctly yields the 
upward convex adiabatic magnetization curve in the or- 
dered region, whereas the mean-field theory yields a flat 



straight line. Another important feature that the spin 
wave theory predicted 17 was the the linear dispersion of 
the magnon in the ordered state. However, the critical 
behavior was not discussed extensively in early studies. 
The HF approximation with the bosonic representation 
predicts that the phase boundary in the T — H plane has 
an algebraic singularity near zero temperature with the 
exponent = 3/2, in contrast to the experimental esti- 
mates. It was reported 53 that by assuming a relativistic 
form for the dispersion relation of magnons the phase 
boundary can be made closer to those estimated exper- 
imentally in an intermediate temperature region, while 
the asymptotic behavior in the zero-temperature limit 
remains the same. A similar observation was made in an 
improved HF calculation 54 based on a realistic dispersion 
relation determined by experiments. 

The BEC mechanism was confirmed by an inelastic 
neutron scattering experiment. 55 A linear dispersion of 
magnons (i.e., tripletons) was observed in the ordered 
phase, as predicted by the spin wave theory 17 and later 
also by bond-operator formalism. 56 This was the first ex- 
perimental observation of the magnon dispersion in the 
ordered phase, while a similar preceding experiment 45 
was carried out in the disordered phase. In the language 
of Boson problems, the linear dispersion can be inter- 
preted as a Goldstone mode characteristic to the BEC. 
This piece of evidence was taken as a clear indication of 
the BEC nature of the phase transition. 

a first-order phase transition at the condensation point 
in contrast to the standard scenario expected for the 
BEC. The reason for this apparent disagreement has not 
been clarified. 

Even without the magnetic field, the BEC phenom- 
ena can be seen when a sufficiently strong pressure 
is imposed. An elastic neutron scattering experiment 57 
demonstrated the existence of the pressure- induced mag- 
netic ordering. The essential mechanism seems to be the 
same as that of the field-induced transition; the compe- 
tition between the gap and the inter-dimer interaction. 
Therefore it is expected that the resulting H — p phase 
diagram at T = is simple and consists of two phases; 
the disordered phase including the origin H = p = and 
the ordered phase surrounding the ordered region. How- 
ever, the phase diagram obtained by the magnetization 
measurement 58 seems to consist of three phases, though 
the very existence of the extra phase or its identity has 
not been established yet. 

While most of the compounds studied in conjunction 
with the field-induced magnetic ordering are quasi one- 
dimensional systems, the compounds BaCuSi206 offers 
an example of a quasi two-dimensional system. An in- 
elastic neutron scattering experiment was carried out 59 
to reveal the two-dimensional nature in the magnon dis- 
persion. Recent heat capacity and magnetization mea- 
surements 60 clearly showed the lambda peak in the heat 
capacity. 

The effect of disorder is certainly an interesting sub- 
ject to explore. It has been already studied through 
the magnetization measurement 50,61 of the compound 
(Tli_ x K x )CuCl 3 . At a finite temperature, the presense 
of the disorder showed only a quantitative effect, namely, 
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reduction of the transition temperature 61 and did not al- 
ter the nature of the phase transition in a qualitative way. 
Qualitative differences could be observed, however, as the 
temperature was lowered. For instance, the asymptotic 
behavior of the critical magnetic field seemed to deviate 
from that of the homogeneous case. Namely, the value 
of the exponent (f> appeared to change when the disorder 
was introduced, as suggested in the scaling theory. 62 

Another type of disorder, namely non-magnetic im- 
purities, was studied for TltCux^Mg^Clg. 63 ' 64 Upon 
Mg doping, a magnetic phase transition is induced even 
at zero magnetic field due to the presense of uncoupled 
magnetic ions. It was found that this impurity-induced 
magnetic ordering does not affect the dispersion of the 
magnons at zero magnetic field. However, experiments 
in the magnetic field are yet to be performed in order 
to observe the effect upon the critical phenomena at the 
QCP. 

4. Scaling and Numerical Results 

As stated above, the bosonic theory explains most of 
the characteristics observed in experiments. However, it 
still has a short-coming, i.e., it is destined to produce 
wrong critical properties for the phase transition at finite 
temperature and there is no systematic way of amend- 
ing the flaw. For example, it results in a discontinuous 
drop in the magnetization at the critical point as we in- 
crease the temperature, whereas we only see a dip in 
experiments. 47 Another apparent disagreement is of the 
estimate of the exponent <j>. Experiments seemed to in- 
dicate a larger value for the exponent <f) for the QCP 
than the HF theory does. Since the incorrectness of the 
mean-field type theory concerning the critical exponent 
is well-known, the disagreement seems trivial at a first 
glance. However, it is not trivial and in fact the HF the- 
ory may be correct concerning the QCP as we see below. 

The first quantum Monte Carlo simulation for clar- 
ifying the critical properties 65 was performed for the 
model Hamiltonian (1) representing weakly coupled lad- 
ders. The estimates of the critical exponent <p differed 
from the HF value; <j> = 1.8(2) for the upper critical field 
and </> = 3.1(2) for the lower critical field. Another Monte 
Carlo simulation was carried out 66 for a similar system. A 
temperature-dependent effective exponent 4>{T) was de- 
fined so that it characterizes the H c — T curve in a finite 
temperature-range centered at T. Again, 4>{T) turned out 
to be greater than 3/2. However, it decreases as the tem- 
perature is lowered, and it was concluded that the trend 
is consistent with the HF value. 

Recently, we carried out another Monte Carlo simu- 
lation 67 to further clarify the problem concerning the 
critical properties at QCP. Since the critical properties 
of the effective Hamiltonian (2) derived with the Tachiki 
and Yamada's mapping is expected to be the same as 
those of the original Hamiltonian (1), we chose to work 
on the effective Hamiltonian for simplicity. The effective 
spin Hamiltonian is a uniform 5 = 1/2 XX Z model in 
three dimensions with a spatial anisotropy and with the 
easy-plane spin anisotropy. Since the spatial anisotropy 
and the S z — S z couplings should not affect the critical 
properties, we can drop them without changing critical 
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properties. Therefore, we arrive at the XY model Hamil- 
tonian. 

n = - jJ2 ( s ? s ! + s ! s j) -hY. s I ( 4 ) 

(if) 1 

The simulation method employed is based on the 
directed-loop algorithm, 68 for which a review can be 
found in an recent article. 69 The magnetization at zero 
temperature was computed as a function of the magnetic 
field. It was found that in the thermodynamic limit the 
deviation of the magnetization from the saturation value 
depends linearly on the deviation of the magnetic field 
from the critical value: 1/2 — (Sf) oc H c \ — H where 
H < H c i. This is consistent with the HF prediction and 
also with the experiments. 

For the susceptibility of the perpendicular components 
of spins, x = J2r la d.T(S x (r, t)S x (0, 0)), the result sug- 
gested that the system size dependence at the critical 
point H = 3 J was characterized by x{fl/L 2 7 H = H c ) oc 
L 2 - 5 . As we see below, this is also consistent with the 
scaling theory that produces the mean-field type scaling 
properties. 

The temperature dependent critical field H C (T) was 
also computed. A set of values of temperature were cho- 
sen and for each temperature the susceptibility was com- 
puted as a function of the magnetic field. We then carried 
out the finite-size scaling analysis with a fixed tempera- 
ture to obtain the critical field at the temperature. The 
finite size scaling plot worked nicely with the critical ex- 
ponents estimated for the classical XY model, namely, 
v = 0.67155(27) and r\ = 0.03 80(4). 70 In this fashion, we 
estimated the critical field at various temperatures. A 
remark may be useful here concerning the critical prop- 
erties at a finite temperature. There has been a confusion 
about the critical behavior of the specific heat, which was 
considered to be divergent in some early articles. In fact, 
the specific heat exponent a = — 0.0146(8) 70 is negative, 
i.e., it has a sharp cusp but finite even at the critical 
point. 

In Fig. 1, two curves are plotted for comparison. The 
solid curve represents the mean- field critical exponent, 
= 1.5, whereas the dashed curve the previous estimate 
<p = 2.0 48 based on an experiment on TICUCI3. At a first 
glance, it seems that the curve with <j> = 2.0 fits the data 
better than <f> = 1.5. However, when the logarithmic scale 
is used, as we do in the inset of Fig. 1, it is rather clear 
that <f) = 2.0 explains only a transient behavior, and the 
correct asymptotic value of the exponent is = 1.5. 

Now, let us see why the mean-field type scaling should 
hold for the QCP in the present system. To this end, we 
first derive the bosonic Hamiltonian from (4) by using 
the mapping Sf — b\, S~[ = bi, and Sf = hi — \ with 
the constraint 

h t = b\bi = 0, 1. (5) 

Then we obtain 

H=-t^2 ( b l b j + b ] b i) -^hi+A^hiihi-l), 

(ij) i i 

(6) 
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Fig. 1. Critical field as a function of temperature. The solid curve 
corresponds to <j> = 1.5 whereas the dashed one <j> = 2.0. The 
inset is the same data plotted in the logarithmic scale. While in 
the intermediate temperature region <f> = 2.0 seem to fit the data 
well, in the low-temperature region, the correct slope <f> = 1.5 
yields a better fitting. 

with new parameters t and /i corresponding to the old 
ones as t ~ J/2 and /i ~ H . The A-term imposes the con- 
dition (5), approximately. However, we believe this ap- 
proximation does not affect the critical properties since 
the difference between the hard (original) constraint and 
the soft one is significant only when more than double 
occupancies occur frequently, for which the average oc- 
cupation number needs to be large. As we approach the 
QCP, however, we have smaller particle density and in 
the end we have zero density at the QCP. Therefore, we 
expect the softening of the constraint has no effect on 
the critical properties. 

The rest of the argument is essentially based on the 
results obtained earlier 62 ' 71 for diluted bose gas. Simi- 
lar results can be found also in textbooks. 72 However, 
the connection to the present model has not been very 
transparent, which we can see in the fact that those early 
scaling theories have been seldom mentioned in conjunc- 
tion with the field- induced ordering in dimer systems. 
Therefore we reproduce the scaling results here, and it 
would be useful to have explicit formulas in the form di- 
rectly comparable with the Monte Carlo results. (In par- 
ticular, explicit formulas for the system-size dependence 
are hardly found in the literature.) First, the continu- 
ous field theory for the model (6) is characterized by the 
action 

S = j dx j dr (V^ + IWI 2 

-/iRei/'-r|V>| 2 + w|V'| 4 ^ ■ (7) 

Obviously, the fixed point is given by h — r = u = 0. A 
simple dimensional analysis suffices to obtain the scaling 
dimensions at this fixed point. The effect of the scale 
transformation around the fixed point up to the scale b 
is the following: 

(3->f3 = (3b- 2 , L^L = Lb' 1 , ip -»• $ = #2, 

h^h = hb 2+ i , r — > f = rb 2 , u — > u = ub 2 ~ d , 



where L is the system size. Note that the non-linearity 
represented by parameter u diminishes as we renor- 
malize for d > 2. Hence, the upper critical dimension 
d c = 2, 71 and the critical properties controlled by the 
u = fixed point when d > 2. For three dimensions, 
therefore, we should expect a mean-field type scaling. 
To be more concrete, let us derive the scaling proper- 
ties for the magnetization, m ~ ~ (Sf). We as- 
sume that a scaling property holds for the singular part 
of $ = — log Z(h, r, u, (3, L), namely, <&(h,r,u, (3, L) = 
$(/i, f, u, (3, L) for an arbitrary b. 

In order to obtain the asymptotic behavior of the crit- 
ical magnetic field, we note that we can make u small by 
choosing large b so that the perturbation theory in u is 
nearly exact for the renormalized system. The perturba- 
tion theory yields that $ for the renormalized system has 
a singularity when the condition 72 f = f c = Au/(3 d ^ 2 is 
satisfied with A being some numerical constant. It means 
that in terms of the bare coupling constants the singu- 
larity occurs when 71 

H c (T)-H c {0)(xT^ 2 . 

This is also confirmed by the simulation. 

5. Summary 

We have presented a brief review of the field-induced 
magnetic ordering. From the early stage of the research, 
the existence of the ordering was established, which was 
also supported by the mean-field theory. The other mean- 
field theory, i.e, the HF theory was later constructed 
based on the bosonic representation and described vari- 
ous features observed in the experiments more precisely. 
Many of recent experimental results, however, seemed to 
contradict even the HF theory in one important point, 
i.e., the critical property at the QCP. There could be 
some possible sources of this discrepancy, such as a pres- 
ence of Dzyaloshinsky-Moria interactions, a staggered 
component of the g-tensor giving rise to an effective stag- 
gered field, or a crystal field that reduces the symme- 
try. However, the recent Monte Carlo simulation results 
present an explanation, though not the unique possibil- 
ity, that the true asymptotic behavior is consistent with 
the scaling theory while the transient behavior is hard to 
get rid of. 

Interesting remaining problems concerning the field- 
induced magnetic ordering are (i) experimental and nu- 
merical investigation of the complete phase diagram of a 
disordered magnet that has a finite gap in the pure case, 
and the critical phenomena at the bose-glass to super- 
fluid transition, and (ii) the effect of a symmetry lowered 
by the other kinds of interactions (crystal fields, effective 
staggered field, etc) or, in case where the gap is caused 
by the easy-plane axial anisotropy, by a transverse field 
tangent to the easy plane. 
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